Capillary-Gravity Waves Generated by a Sudden Object Motion 
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We study theoretically the capillary-gravity waves created at the water-air interface by a small 
object during a sudden accelerated or decelerated rectilinear motion. We analyze the wave resis- 
tance corresponding to the transient wave pattern and show that it is nonzero even if the involved 
velocity (the final one in the accelerated case, the initial one in the decelerated case) is smaller 
than the minimum phase velocity Cmin = 23 cms -1 . These results might be important for a better 
understanding of the propulsion of water- walking insects where accelerated and decelerated motions 
frequently occur. 

PACS numbers: 47.35.-i,68.03.-g 



I. INTRODUCTION 

If a body (like a boat or an insect), or an external 
pressure source, moves at the free liquid-air interface, it 
generates capillary- gravity waves. These are driven by a 
balance between the liquid inertia and its tendancy, un- 
der the action of gravity and under surface tension forces, 
to return to a state of stable equilibrium [l[ . For an in- 
viscid liquid of infinite depth, the dispersion relation of 
capillary-gravity waves, relating the angular frequency w 
to the wavenumber k is given by uj 2 — gk + ^k 3 / p, where 
7 is the liquid-air surface tension, p the liquid density, 
and g the acceleration due to gravity @. The energy 
carried away by the waves is felt by the body (or the 
pressure source) as a drag K w , called the wave resistance 
Q. In the case of boats or ships, the wave resistance 
has been well studied in order to design hulls minimizing 
it The case of objects that are small compared to 
the capillary length n~ 1 — y^/Jfig) has been considered 
only recently [a-Hlj. 

In the case of a disturbance moving at constant veloc- 
ity V = |V| on a rectilinear trajectory, the wave resis- 
tance K w is zero for V < c m in where c m i n — (^gp/^f) 1 ^ 4 
is the minimum of the wave velocity c(fc) = uj{k)/k — 
y 'g/k + jk/p for capillarity gravity waves 0, M, 03 • Ef- 
fectively, in the frame moving with that object, the prob- 
lem must be stationary. That is true only if the ra- 
diated waves have a phase velocity c(k) equal to the 
object's velocity V. If V < c min , no solutions exist, 
hence no waves and the wave resistance is zero. For wa- 
ter with 7 = 73mNm -1 and p = 10 3 kgm -3 , one has 
Cmin — 0.23 ms -1 (at room temperature). It was re- 
cently shown by Chepelianskii et al. [l3| that no such 
velocity threshold exists for a steady circular motion, for 
which, even for small velocities, a finite wave drag is ex- 
perienced by the object. Here we consider the case of a 
sudden accelerated or decelerated rectilinear motion and 
show that the transient wave pattern leads to a nonzero 
wave resistance even if the involved velocity (the final 
one for the accelerated case, the initial one for the decel- 



erated case) is smaller than the minimum phase velocity 
Cmin = 23 cms -1 . The physical origin of these results 
is similar to the Cherenkov radiation emitted by acceler- 
ated (or decelerated) charged particles [13, EH ■ 



II. EQUATIONS OF MOTION 

We consider an inviscid, deep liquid with an infinitely 
extending free surface. To locate a point on the free sur- 
face, we introduce a vector r = (x,y) in the horizontal 
plane associated with the equilibrium state of a fiat sur- 
face. The motion of the disturbance in this plane induces 
a vertical displacement £(r, t) (Monge representation) of 
the free surface from its equilibrium position. 

Assuming that the liquid equations of motion can be 
linearized (in the limit of small wave amplitudes), one 
has fil 



d 2 C(Kt) 
dt 2 



w (fc) 2 C(M) 



p 



(i) 



where P ext (k, t) and C(k, t) are the Fourier transforms of 
the pressure distribution and the displacement, respec- 
tively 16(. In what follows, we will assume that the pres- 
sure distribution is axisymmetric around the point r (t) 
(corresponding to the disturbance trajectory). P ex t(k, t) 
can then be written as P ex t(k)e~ l]s " rotyt \ 



A. Uniform straight motion 

Let us first recall the results previously obtained in the 
case of a uniform straight motion 0, 0, H3, 13 • Such a 
motion corresponds to r (t) = (— Vt, 0) where V is the 
constant velocity of the disturbance. 

Equation (flj then becomes 
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+ L0(k) 2 ((k,t) 
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FIG. 1: Wave resistance R™ (in units of Poft/7) as a func- 
tion of the reduced velocity V/cmin for a uniform straight 
motion, see Eq.([5]). The pressure disturbance is assumed to 
be a Lorentzian, with a Fourier transform P ex t(k) = Poe~ bk , 
where b is the object size (set to b — 0.1k _1 ). 
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FIG. 2: The wave resistance R w (in units of Pq Jt/7) is shown 
as a function of the reduced time Cminut for an accelerated 
motion with different U = V/c m in (see Eq.©). Respectively, 
panels (a), (b), (c) and (d) correspond to a reduced velocity 
U = 1.5, 1.1, 0.9 and 0.5. 



The above equation corresponds to the equation of an 
harmonic oscillator forced at angular frequency Vk x . 
We can solve it by looking for solutions with a time- 
dependence of the form e tVk ^ t ^ leading to 



C(M) 



kP ext (k) 



JVkxt 



P (w{ky - {vk x f) 



(3) 



Following Havelock Q, the wave resistance R w experi- 
enced by the moving disturbance is given by 



R„ 



f f dk x dk y 



kC(M) P*(M). 



(4) 



This expression for the wave resistance represents as 
the total force exerted by the external pressure on the 
free surface R w = JJ dxdyP ext {r,t)S7((r,t) written in 
Fourier space. Using Eq.Q and integrating over the an- 
gular variable one obtains [5j 



B. Accelerated straight motion 

We now turn to the case where the disturbance - 
initially at rest - is suddenly set to a uniform mo- 
tion (characterized by a constant velocity V) at time 
t = 0. The corresponding trajectory is given by i"o(t) = 
— V 1 6(t) u x . As long as the perturbation does not move 
(i.e. for t < 0), the wave resistance is equal to zero. In 
order to calculate the wave resistance for t > 0, we solve 
Eq.([2]) along with the initial conditions £(k, t = 0) = 

— = 0, yielding 



and 



C(k, t) = - fdT kP ^ k) e-^) sin (w {k) {t _ T) y 
Jo u{k)p 

(6) 

Equation (0} then leads to the following expression for 
the wave resistance 



R„ 



D dkkP[ 
2^ ~V 



_ 2bk 0(V-c(k)) 
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(5) 



where 9(x) is the Heaviside step function and u x the 
unit vector along the x-axis. The behavior of as a 
function of the disturbance velocity V is illustrated in 
Fig. [TJ There we have assumed a pressure disturbance 
of Lorentzian form with a Fourier transform P ext (k) = 
Poe~ bk , where b is the object size (set to 6 = O.Ik -1 ). 
This choice will be taken for the figures throughout this 
work. The wave resistance is equal to zero for V < c m i n 
and presents a discontinuous behavior at V — c m i n (see 
reference @ and [l3[ for a more complete discussion). 



R w (t) 
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where J\{x) is the first Bessel functions of the first kind. 
In the long time limit, one has [19( 



lim / sin(w (k) u) J\(kVu) du 



uj{k)9{V - c(k)) 



V 2 k 2 y/1- (c(k)/V) 2 
(8) 

and therefore the wave resistance Eq. ([7]) converges to the 
uniform straight motion result, Eq. ([5]). The behavior of 
Ru,(£) is represented on FigJ5]for different values of the 
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disturbance velocity and will be discussed in detail in 
section IIIIl 



C. Decelerated straight motion 

Let us now consider the case of a disturbance moving 
with a constant velocity V for t < 0, and suddenly set 
and maintained at rest for t > 0. This corresponds to the 
following trajectory: r (t) = — V 1 6(— t) u x . As long as 
t < 0, the wave resistance R m is given by the uniform 
straight motion result Eq. (JSJ) - In order to calculate the 
wave resistance for t > 0, we solve Eq.© along with the 
initial conditions 



C(M = o) 



and 



at 



(k,t = 0) = 



p(w(A;) a - (W*) 2 ) 
p(w(fc) 2 - (Ffc^) 2 ) 



(9) 



(10) 



(where we have used Eq. JSJl). This leads to 



C(M>o) = 

( kP ex t{k) k P ex t(k) 



cos (w (fe) i) 



^ p W (fc) 2 p(w(fc) a - (W*) 2 ), 

k (iVk x ) P ext (k) sin (a; (k) t) _ kP ext (k) 
p{uj{kf - {Vk x ) 2 ) u(k) poj(k) 2 ' 



(11) 



Equation ((4|) then leads to the following expression for 
the wave resistance: 

Ruj(i) = 

dk k\P ext (k)f 9{V-c{k)) 
p.V- C0S M fc )*) , 



D dk k\P ext (k)\ 2 
o 2tt p.V 2 



sin (oj (k) t) —— — 




In the l ong time limit (t — ¥ oo), the Riemann-Lebesgue 
lemma [201 ] . for a Lebesgue integrable function / 



lim / f(x)e lxt dx 

t— >oo 







(13) 



permits to determine the limit: the wave resistance given 
by Eq. (fT2"|) converges to for t — > oo, as expected. The 
behavior of Ru,(f) is represented in FigJ3] for different 
values of the disturbance velocity and will be discussed 
in detail in section the next section (Sec JIIIj) . 
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FIG. 3: The wave resistance K w (in units of Pq k/j) is shown 
as a function of the reduced time c m i„Kt for a decelerated 
motion with different reduced velocities U = V/c m in (see 
EQ. (fr2|) '). Respectively, panels (a), (b), (c) and (d) correspond 
to a reduced velocity U = 1.4, 1.1, 0.98 and 0.7. 



III. RESULTS AND DISCUSSION 

A. Accelerated straight motion 

Figures [5] and [3] represent the behavior of the wave 
resistance for the accelerated and the decelerated cases, 
respectively. In order to get a better understanding of 
the behavior of R„, = • u x , we will perform ana- 
lytic expansions of Eq.© and Eq. flT2")) . respectively. Let 
us start with the accelerated case, Eq.©. Since one 
has the product of two oscillating functions, sin (w (k) u) 
and Ji (kVu), one can use a stationary phase approxi- 
mation [2l|. The sine function oscillates with a phase 
cf>i = Lo(k)u whereas the Bessel function J\ oscillates 
with a phase 4>2 — Vku. Their product has thus two 
oscillating terms, one with a phase (f>— = 4>\ — 4>2 and 
the other with a phase 0+ = 4>\ + 4>2- According 
to the stationary phase approximation, the important 



wavenumbers are given by 



and 



0. 



dk dk 
The latter equation does not admit any real solution 

and the corresponding contribution to the wave resis- 
tance decreases exponentially and can be neglected. The 



equation 

Cg(k) = 



dk 



= leads to {c g (k) — V)u — 0, where 

d ^ fc ^ is the group velocity of the capillary- 
gravity waves [3J. This equation has solutions only if 
V > min(c ff (fe)) = (3^3/2 - 9/4) 1 / 4 c mm w 0.77c mm . 
When this condition is satisfied, two wavenumbers are se- 
lected: kg (mainly dominated by gravity) and k c (mainly 
dominated by capillary forces), with k g < k c (see FigU]). 
If these two wavenumbers are sufficiently separated (that 
is, for velocities not too close to 0.77c m i„), one then finds 
that in the long time limit R^ oscillates around its final 
value as 
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FIG. 4: Graphical representation of wavenumbers fci, fc2, k c 
and k g in units of k. k\ and k^ are the solutions of the equa- 
tion c(k) = V and correspond to the intersection between the 
curve c{k)/c m in and the line U = V/c min . k g and k c are 
the solutions of the equation c g (k) = V and correspond to 
the intersection between the curve c g (k) /c m in and the line 
U = V/cmin- Analytical expresions of k c and k g can be given 
but are rather lengthy. 



FIG. 5: The wave resistance R w (in units of Pq Jt/7) is shown 
as a function of the reduced time Cminut for an accelerated 
motion with different U = V/c m in (see Eq.©). Respectively, 
panels (a), (b) and (c) correspond to a reduced velocity U = 
1.5, 1.1 and 0.9. The solid lines are obtained by a numerical 
integration of Eq.([7|. The dashed lines correspond to the 
asymptotic expansion given by Eq. (|14p . 



R w (t) = R^(oo) 



1 



kc /2 \P ex t(k c )\ 2 cos(fi c t) 



-•(/-. >\/||£(M 



2tt P VV 

1 k 5 g /2 \P ext (k g )\ 2 sm(n g t) 



flrt 



2tt P VV 



(14) 



where R tu (oo) is given by Eq.© (and is equal to zero if 
V < Cmin), Q c = (c(fc c ) - c g (k c ))k c and fl g = {c(k g ) - 
c g(kg))k g - 

Therefore, even if the disturbance velocity V is smaller 
than c m i n , there exists a transient nonzero wave resis- 
tance decreasing as 1/t (for V > 0.77 c m in)- 

We obtain a good agreement between the numerical 
calculation of Eq.© and the analytical approximation 
Eq. (fT4")) as shown on FigJS] Note that the oscillations 
displayed by the wave resistance are characterized by a 
period 2n/Q, c or 2ir/il g that diverges as V approaches 
c m in- In the particular case where V ^> c m in, Eq. fT4"|) 
reduces to 



p AiA T> f n „\ i 1 KC min\Pe-xt(k g )\ 2 fc 

R w (t) - R u) (oo) + ^ 75 - — sm^- 

2*/ 2 K\P ext {k c )\ 2 f8V 3 Kt 



„Kt 



jpVt 



-cos 



27c? 



8V 



(15) 



Let us now give a physical interpretation for the wave 
resistance behavior as described by Eq. ([T4")) : during the 



sudden acceleration of the disturbance (taking place at 
t = 0), a large range of wavenumbers is emitted. Waves 
with wavenumbers k such that c g {k) > V will move faster 
than the disturbance (which moves with the velocity V) , 
whereas waves with wavenumbers k such that c g (k) < V 
will move slower. The main interaction with the mov- 
ing disturbance will therefore correspond to wavenum- 
bers satisfying c g (k) — V, that is to k c and k g (hence 
their appearance in Ea.(fI3])). Due to the Doppler effect, 
the wave resistance R^, (which is the force exerted by the 
fluid on the moving disturbance) oscillates with an an- 
gular frequency (c(fc) — V)k, hence the appearance of Q c 
and ft g in Eq.(IT4l. 

Note that in the case V < 0.77c m in, the wave re- 
sistance, Eq.©, is nonzero and decreases exponentially 
with time (see Fig|2{d)). 

In order to get a better physical picture of the gen- 
erated wave patterns, we have also calculated numeri- 
cally the transient vertical displacement of the free sur- 
face C( r >0 m the accelerated case (Eq.©). The corre- 
sponding patterns (as seen in the frame of the moving 
object), are presented on Figs. El [7] and [8] for different 
reduced times c m inKt (1, 10 and 50, respectively). At 
CminKt = 1, the perturbation of the free surface is very 
localized around the disturbance (close to the same one 
obtained by a stone's throw). At c m i n Kt = 10, some 
capillary waves can already be observed at the front of 
the disturbance. At c m in,Kt — 50, on can see a V-shaped 
pattern that prefigure the steady pattern of the uniform 
straight motion (obtained at very long times). These pre- 
dictions concerning the wave pattern might be compared 
to experimental data using the recent technique of Moisy, 
Rabaud, and Salsac [23j. 
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FIG. 6: Accelerated straight motion: Transient vertical dis- 
placement (in units of Po k / (pc^in)) of the free surface at 
tc m i n n = 1 obtained by inverse Fourier transform of Eq.©. 
Note that the surface disturbance is localized around the ob- 
ject and close to the one obtained by a stone thrown in water. 
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FIG. 7: Accelerated straight motion: Transient vertical dis- 
placement (in units of PoK/(pc^ n i n )) of the free surface at 
tCminK, — 10 obtained by inverse Fourier transform of Eq.((6|. 
One can already see ahead of the disturbance the waves as- 
sociated with fci, while the ones associated to fe are less well 
formed. 
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FIG. 8: Accelerated straight motion: Transient vertical dis- 
placement (in units of PoK/(pc^ n i n )) of the free surface at 
tcminn = 50 obtained by inverse Fourier transform of Eq.((6|. 



B. Decelerated straight motion 

Let us now turn to the case of the decelerated mo- 
tion, described by Eq. (fl~2|) . We first will assume that 
V > c m i n . In that case, the denominators appearing in 
the integrals vanish for k such that c(k) = V, that is for 
k\ = k(U 2 — \ZU 4 — 1) (mainly dominated by gravity) 
and = n(U 2 + s/U 4 — 1)) (mainly dominated by capil- 
lary forces) as shown in FigfJ] where U — V/c min . Both 
wavenumbers contribute to the wave resistance. Using 
the stationary phase approximation one then finds that in 
the long time limit R w oscillates as cos {Vk\t + 7r/4) / \/t. 
More precisely, one has (for large t) 



R w (t) = ' — cosCW + 7r/4). 

(16) 

The wave resistance displays oscillations that are char- 
acterized by a period 2ir/ (Vki). Figure H] shows the good 
agreement between the numerical calculation of Eq. (fT2|) 
and the analytical approximation Eq.(|16j) for large times 
{c m i n K,t > 10). Again one can give a simple physical 
interpretation of the wave resistance, Eq. pl))) . The pe- 
riod of the oscillations is given by 27r/(Ufci) and depends 
only on the wavenumber ki (mainly dominated by grav- 
ity), even in the case of an object with a size b much 
smaller than the capillary length. This can be under- 
stood as follows. For t < 0, the disturbance moves with 
a constant velocity V and emits waves with wavenumber 
k\ and waves with wavenumber &2 (mainly dominated 
by capillary forces). The fci-waves lag behind the distur- 
bance, while the fc2-waves move ahead of it. When the 
disturbance stops at t = 0, the /c 2 -waves keep moving for- 
ward and do not interact with the disturbance. However, 
the fci-waves will encounter the disturbance and interact 
with it. Hence the period 2ir/(yk\) of the wave resis- 
tance Eq. ([T5| . The \j\fi decrease of the magnitude of 
wave resistance in Eq. (fTl)|) can be understood as follows. 
At time t > 0, the disturbance (which is at rest at x = 0) 
is hit by a /ci-wave that has previously been emitted dis- 
tance c g (ki)t away from it. The vertical amplitude C of 
this wave is inversely proportional to the square root of 
this distance. Indeed, as the liquid is inviscid, the energy 
has to be conserved and one thus has ( cx 1/ y / c g (ki)t. 
Since the wave resistance R^, is proportional to vertical 
amplitude of the wave £, on recovers that H w decreases 
with time as 1/ yjc g {ki)t. 

In the particular case where V ^> c min , Eq. dT^l) re- 
duces to: 



R w (t) 



1 K 3 / 2 \P ext (k c )\ 



2 3 
rain 



2^ P V n / 2 Vi 



cos 



-t + 7T/4 



2V 

Let us now discuss the case V < c m i n (still considering 
the decelerated motion Eq . (TEH) ) . In that case, the wave 
resistance oscillates (in the long time limit) approxima- 
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FIG. 9: The wave resistance K w (in units of Pqk/j) is shown 
as a function of the reduced time CminKt for a decelerated 
motion with different reduced velocities U = V/c m in (see 
Eg. ()12|l V Respectively, panels (a), (b), (c) and (d) corre- 
spond to a reduced velocity U = 1.5, 1.01, 0.98 and 0.7. The 
dashed lines correspond to the asymptotic expansion given by 
Eq.CEHl). 



tivelyas (e UV1 UHKlCmin / yjtnc mil ^j sin(U 3 tKc mm + \)- 
More precisely, one has (for large t) 



where 



2 K 2 \P e xt (k)\ 2 sin([/ 3 c mm nt + x) 



-UVl-U A tKc min 



V27rC/V2(i-[/4)i/4( 3f/ 4 + 1) i/4' ( 18 ) Sa l sac t 



X = (3/2) arctan Ufl - U A /U 2 ^j 

- (1/2) arctan Ufl - U 4 /(2U 2 )^j (19) 



and k = k (U 2 + i\J\ — U 4 ) . Note that in the case 

V < c m i n , the wave resistance Eq. (fl~B|) also displays some 
oscillations, although no waves were emitted at t < 0. 
These oscillations - that have a exponential decay in 
time - might be due to the sudden arrest of the distur- 
bance. Such oscillations could also be present in the case 

V > Cmin but are hidden by the interaction between the 



fci-waves and the disturbance which lead to a much more 
slower l/^/t decay. 



IV. CONCLUSIONS 

In this article, we have shown that a disturbance un- 
dergoing a rectilinear accelerated or decelerated motion 
at a liquid-air interface emits waves even if its velocity 
V (the final one in the accelerated case and the initial 
one in the decelerated case, respectively) is smaller than 
Cmin- This corroborates the results of Ref.[l3j. For this 
purpose, we treat the wave emission problem by a lin- 
earized theory in a Monge representation. Then, we de- 
rive the analytical expression of the wave resistance and 
solve it by numerical integration. Asymptotic expansions 
permit to extract the predominant behavior of the wave 
resistance. Some vertical displacement patterns are also 
calculated in order to shown how the waves invade the 
free surface. 

The results presented in this paper should be im- 
portant for a better understanding of the propulsion 
of water-walking insects like whirligig beetles, 

where accelerated and decelerated motions frequently oc- 
curs (e.g., when hunting a prey or escaping a predator 
|28j). Even in the case where the insect motion appears 
as rectilinear and uniform, one has to keep in mind that 
the rapid leg strokes are accelerated and might produce a 
wave drag even below c m i n . The predictions concerning 
the wave patterns might be compared to experimental 
data u sing the recent technique of Moisy, Rabaud, and 
Salsac [231 ] . 

It will be interesting to take in our model some non- 
linear effects [2!| because the waves radiated by whirligig 
beetles [28[ have a large amplitude. Recently Chepelian- 
skii et al. derived a self-consistent integral equation de- 
scribing the flow velocity around the moving disturbance 
[30| . It would be interesting to incorporate this approach 
into the present study. 
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